FACTORIZATION OF THE HYPERGEOMETRIC-TYPE 
DIFFERENCE EQUATION ON THE UNIFORM LATTICE 

R. ALVAREZ-NODARSE*, N. M. ATAKISHIYEVt , AND R. S. COSTAS-SANTOS* 

Abstract. We discuss factorization of the hypergeometric-type difference equations on the 
uniform lattices and show how one can construct a dynamical algebra, which corresponds to each of 
these equations. Some examples are exhibited, in particular, we show that several models of discrete 
harmonic oscillators, previously considered in a number of publications, can be treated in a unified 
form. 



1. Introduction. The study of discrete system has attracted the attention of 
many authors in the last years. Of special interest are the discrete analogs of the 
quantum harmonic oscillators [l[51IHl[Tni[TIl[Tni[T71[THl[lTl[25] among others. 

There are several methods for studying such systems. One of them is the fac- 
torization method (FM), first introduced for solving differential equations [TH]. This 
classical FM is based on the existence of the so-called raising and lowering operators 
for the corresponding equation, which allow to find the explicit solutions in a simple 
way, see e.g. [6l[22]. Later on, Miller extended it to difference equations (23j and 
g-differences -in the Hahn sense- [23] . In the case of difference equations this method 
has been also extensively used during the last years (see e.g. [SI UHl [131 [HI HH] for 
difference analogs on the uniform lattice and [21 [31 [3 [El [IHl [HI [13] for the q-case) . 

Later on, references [HI 13 [H] indicated a way of constructing the so-called "dy- 
namical symmetry algebra" by applying the FM to differential or difference equations 
[H [ini [TT] and then this technique has been used to consider some particular in- 
stances of g-hypergeometric difference equations. Of special interest is also the paper 
by Smirnov [55], in which the equivalence of the FM and the Nikiforov et al formu- 
lation of theory of g-orthogonal polynomials [13], was established. In [5], following 
the papers [T31 [21] for the classical case, it has been shown that one can factorize the 
hypergeometric-type difference equation (12. ip in terms of the above-mentioned raising 
and lowering operators. 

Our main purpose here is to show how to deal with all different cases of difference 
equations on the uniform lattice x{s) = s in an unified form. One should consider 
this paper as an attempt to provide a background for the more general q-linear case 
(since in the limit as q goes to 1, the g- linear case reduces to the uniform one). Some 
results concerning this general case will be also given in the last section. 

The structure of the paper is as follows. In Section 2 some necessary results on 
classical polynomials are collected. In section [3] the factorization of the hypergeome- 
tric-type difference equation is discussed, which is used in section [3] to construct a 
dynamical symmetry algebra in the case of the Charlier polynomials. In section [5] 
the Kravchuk and the Meixner cases are considered in detail. Finally, in section [5] we 
briefly discuss a possibility of applying this technique to the g-case. 
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FACTORIZATION OF DIFFERENCE EQUATIONS 



2. Preliminaries: the classical "discrete" polynomials. The discretization 
of the hypergeometric differential equation on the lattice x{s) [551 leads to the 
second order difference equation of the hypergeometric type 

, , A yy{x{s)) , ^ , Ay{x{s)) ,,,,,, „ 



Ax(s-l/2) Vx{s) Ax{s) «v 

where A/(s) /(s + 1) - /(s), V/(.s) := f{s) - f{s - 1). 

The most simple lattice is the uniform one x{s) — s and it corresponds to the 
equation 

a{x)A\'y{x) + T{x)Ay{x) + Xy{x) = 0. (2.2) 

The above equation have polynomial solutions P„(s), usually called classical discrete 
orthogonal polynomials, if and only if A = A„ = —n{T' + (n — 1)ct"/2). 

It is well known |26j that under certain conditions the polynomial solutions of 



(|2.2p are orthogonal. For example, if <7{s)p{s)x'' 
the polynomial solutions Pn{x) of (|2.2p satisfy 

b-l 

(P„, Pyn)d = ^ P-n{s\ q) Pmis; q) p{s) 



0, for all fc = 0, 1,2, . . ., then 



(2.3) 



where the weight functions p{s) are solutions of the Pearson-type equation 

A[a{s)p{s)]^T{s)p{s) or a{s + l)p{s + I) ^ [(7{s) + t{s)]p{s). (2.4) 
In the following we will consider the monic polynomials, i.e., Pn{x) — x"+bnx"^^+- ■ ■ . 



The polynomial solutions of ()2.2p are the classical discrete orthogonal polynomials 
of Hahn, Meixner, Kravchuk and Charlier and their principal data are given in the 
table mH 

They can be expressed in terms of the generalized hypergeometric function pFq, 

(ai)fc(a2)fc • • • (ap)fc x'' 



pr q 



fll, 02, . . . , Op 
61,62, • • ■ ,6g 



^ (6i)fc(62)fc---(6g)fe A:! ' 



where {a)k is the Pochhammer symbol (or shifted factorial) 

(a)o = 1, {a)k = a{a + l)(a + 2) • • • (a + fc - 1), fc = 1, 2, 3, . . 



(2.5) 



(2.6) 



Using the above notations, we have for the monic polynomials of Hahn, Meixner, 
Kravchuk and Charlier, respectively 



[a + p + n + l)n 
(7)«M" 



-x,a + (3 + n + 1, 
1 - TV, /3 + 1 



M2'^ix) 
Klix 



{p - 1)« 

{-pYn\ 

' [N -n)\ 



2F1 
2F1 



-n, —X 
7 

— n, —X 
~N 



1 

1 



(2.7) 
(2.8) 
(2.9) 



C^(x) = (-/.)" 2F0 



(2.10) 



A further information on orthogonal polynomials on the uniform lattice can be 
found in m [201 [23 [21] ■ 
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Table 2.1 

The classical discrete orthogonal monic polynomials. 





Hahn 

K'^{x;N) 


Meixner 
M;]:'''(a::) 


Kravchuk 
Kl(x) 


Charlier 

C^(a::) 


[a,h] 


10,N] 


[0,oo) 


[0, iV + 1] 


[0,00) 


g{x) 


x{N + a ~ x) 


X 


X' 


X 


t(x) 


(p + iJ(,Jv — ij — (a + /3 + l)x 




iVp — X 

i-p 


fl — X 


a + T 


(a: + /3 + l)(iV - 


/ix + 7^ 






An 


n(n + a + /J + 1) 


(1 - /i)n 


n 


n 


p{x) 


r(iV+c«-a;)r{/3+a; + l) 


M==r(7+a;) 


P G (0, l),n <Af-l 




r{Ar_a;)r(a; + l) 

a,/3 >-l, n < N-\ 


r(7)r(a:+l) 

7> 0,^ G (0, 1) 


r(a: + l) 

/t > 




n!r(a + /3 + JV + n+l) /a + ^ + 2n + 2\-l 
(]V-n-l)!(c« + ;9 + n + l)„ \ a + n+1 j 


(i-m)-'+^" 







3. Factorization of the difference equation. Let us consider the following 
second order linear difference operator 

[)i(.s) = -iy{s - 1) e-^= - iy{s) + [2a(s) + t(s)]/, (3.1) 

where e"^"7(s) = f{s + a) for aU a e C, i'(s) = ^'^(s + l)[cr(s) + t(s)], and / is the 
identity operator, and let be the set of functions 



$«(s) = ^^P„(s), (3.2) 

where dn is a norm of the polynomials Pn{s), which satisfy equation (j2.2p . and p{s) is 
the solution of the Pearson- type equation (|2.4I) . If P„(s) possess the discrete orthog- 
onality property (j2.3p . then the functions ^n{s) have the property 

6-1 

($n(s), $m(s))d = ^ $„(s)$,„(s) = (5„,™. (3.3) 

s— a 

Using the identity V = A — VA and the equation p.2p . one finds that 

f)i(s)$„(s) = A„$„(s), (3.4) 

i.e., the functions $„(s), defined in (13.21) . are the eigenfunctions of t)i(s). In the 
following we will refer to f)]^(s) as the hamiltonian. 

Our first step is to find two operators a(s) and b{s) such that the Hamiltonian 
i)i{s) — b{s)a{s), i.e., the operators a{s) and b{s) factorize the Hamiltonian i)i{s). 

Definition 3.1. Let a be a real number. We define a family of a- down and a -up 
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operators by 

aiis):^ (e^^V^- V'^i^) + r{s) l) , 

(3.5) 

alis):^ (v^e-^= - Va(s)+r(s)j)e"^% 

respectively. 

A straightforward calculation (by using the simple identity e^-' V = A) shows 
that for all a € M 

f)i(s) = aUs)aiis), 

i.e., the operators Cl^(s) and Cl^(s) factorize the Hamiltonian, defined in p.ip . Thus, 
we have the following 

Theorem 3.2. Given a Hamiltonian [)^{s), defined by \3.1\) . the operators 
Cl^(s) and ClJ,(s), defined in 113.5]} . are such that for all a e C, the relation i)i{s) = 
<^l{s)<^i{s) holds. 

4. The dynamical algebra: The Charlier case. Our next step is to find 
a dynamical symmetry algebra, associated with the operator i)-^(s), or, equivalently, 
with the corresponding family of polynomials, i.e.. To find two operators a{s) and 
b{s), that factorize the hamiltonian f)i(s), i.e., f)i(s) = b{s)a{s), and are such that 
its commutator [a{s),b{s)] = a(s)b{s) — b{s)a{s) = I, where I denotes the identity 
operator. 

Theorem 4.1. Let i)i{s) be the hamiltonian, defined in 113.1]) . The operators 
b{s) — ClJ,(s) and a{s) — Cl4(s), given in 113.5]) . factorize the Hamiltonian f)]^(s) iS. 1]) 
and satisfy the commutation relation [a{s), b{s)] — A for a certain complex number A, 
if and only if the following two conditions hold: 

a{s - a)[a{s - a) + t{s ~ a)] ^ 

a{s)[a{s - 1) + t{s - 1)] ^ ' ' 

and 

a{s - a + 1) + a{s - a) + t(s - a) - 2a{s) - t{s) = A. (4.2) 

Proof. Taking the expression for the operators a^(s) and a^{s), a straightforward 
calculation shows that a^is)a'l,{s) — Ai{s)e^'' + ^2(5)6^'^= + ^3(5)/, where 

Ai{s) = -^/a{s + 1 - a)[a{s - a + 1) + t{s - a + 1), 

^2(5) = -\/cris - a)[a{s - a) + r(s - a), (4.3) 

^3(5) ~ a(s + 1 — a) + (t(s — a) + r(s — a). 
In the same way, al^{s)aj^{s) — [)j(s) = i3i(s)e'^= + B2{s)e^^' + i?3(s)/, where 

Bi(s) = -Ks), -82(5) = -1^(5-1), Bsis)^2a{s)+T{s). (4.4) 
Consequently, 

= {A^{s)-B,{s))e^^ + (^A,{s)~B,{s))e-^^ + [A,{s)-B^{s))L (4.5) 
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To eliminate the two terms in the right-hand side of (j4.5p . which are proportional 
to exp(±9s), one have to require that Ai{s) — i?i(s) = and ^2(3) — -62(3) — 0. 
But Ai{s) / Bi{s) =^2(5 + 1)/B2{s + 1), hence, the requirement that Ai{s) = -Bi(s) 
entails the relation ^2(5) = ^2(3), and vice versa. Thus, from (|4.5p it follows that 
the commutator [al{s) , aj;^{s)] = A, iff Ai{x) = i?i(s) and ^3(5) - Sals) = A. □ 

Using the main data for the discrete polynomials (see table 12. ip , we see that the 
only possible solution of the problem 1 corresponds to the case when a{s) + r(s) = 
const, and a — Oi.e., the Charlier polynomials. Moreover, in this case A„ = n. 

Corollary 4.2. For the hamiltonian, associated with the Charlier polynomials, 

()f (s) = -^/i^e-^'■ - v/(s + l)Me^^ + + t^)I, 



f)f (s)$^(s) = $g(s) ^ y "'217" C^{s). A'>0, n = 0,l,2,.... 

Furthermore, the operators 

a^is) = V^e^'^ -^I, alis) = yje-^= -y^I, (4.6) 

are such that f)f = oJ(s)aQ(s) and [aQ(s), aj(s)] = 1. 
Notice that, since P);^(s)$(s) = A$(s), 

f)i(s){<(s)$(s)} = alis)aiis){atisMs)} = (ai(s)aS(s) - i){ai(s)$(s)} 

^{X-i){at{sMs)}, 
l),{s){al{sMs)} = al{s)at{s)aUsMs) = al{s){X + i)$(s) 

= {X + i){al{sMs)}. 

In other words, if $(s) is an eigenvector of the hamiltonian t)i{s), then aQ(s)$(s) is 
the eigenvector of f)i(s), associated with the eigenvalue A — 1, and aj(s)$(s) is the 
eigenvector of l)i(s), associated with the eigenvalue A + l. In general then [Oq]''(s)$(s) 
and [aj]'"'(s)<l'(s) are also eigenvectors corresponding to the eigenvalues X—k and A+fc, 
respectively. 

Using the preceding formulas for the Charlier polynomials, one finds 

aj(s)$^(s) = ;7„$^+i(s), ai(s)<i>^(s) - A. $^-1(5), (4.7) 

where [/„ and Dn are some constants. 

If we now apply ClJ(s) to the first equation of (|4.7p and then use the second one and 
(|3.4I) . we find that A„ — _D„L/„_i. On the other hand, applying CIq(s) to the second 
equation in (14. 7p and using the first one, as well as the fact that aQ(s)aJ(s)$$^(s) = 
(A„ + l)$^(s), one obtains that 1 + A„ = J7„-Dn+i = A„+i, from which it follows that 
A„ should be a linear function of n (that is also obvious from table [^TT|) . 

If we use the boundary conditions a{s)p{s)\^_^ ^ = 0, as well as the formula of 
summation by parts, we obtain 

(a;^(s)$„(s),$„(s))d = ($™(s),aj(s)$„(s))d. 
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i.e., the operators CIq(s) and ClJ(s) are mutually adjoint. 

From the above equality (the adjointness property) and (|4.7I) it follows that 
Dn+i = Un, thus C/^ = A„+i, therefore Un = \/A„+i and Ai = i-S-, we 

have the following 

Corollary 4.3. The operators ClJ(s) and Clg(s) are mutually adjoint with respect 
to the inner product (•,•)£; and 



From the above corollary one can deduce that 



(4.8) 



v/^<(s + l)- Va^<I>C^(.s) =0 ci^C(^s) = No^^ 

Using the orthonormality of $o'(s), one obtains that A^o — e^^l"^ . Thus 

<i>^(^) - ^ K(.)]"<i>o^(.) = ^ [yJe^= - VaI/]" (/^) ■ 
Notice that 

[t)i(s),aj(s)] - V7^(M-i) + /ia5(s), [fii(s),ai(s)] = -V7^(A*-i)-A*ai(s). (4.9) 

This example constitute a discrete analog of the quantum harmonic oscillator [5] . 

5. The dynamical algebra: The Meixner and Kravchuk cases. From the 
previous results we see that only the Charlier polynomials (functions) have a closed 
simple oscillator algebra. What to do in the other cases? To answer to this question, 
we can use the following operators: 



1 , 



a{s) = ^(t(s + 1) e2<^= - y^a{s-l) + T{s~l)e-'2'^% 

(5.1) 

a+(s)= e^2^=^cr(s + l) - e2^-^cr(s-l) + r(.s-l). 
For this operators 

[)i(s) =a(s)a+(s) + r'-(7". 

We will define a new hamiltonian i)2is) and operators 6(s) and b^{s) 

l)2(s) = C2f)i(s) + A 6(s)-aa(s) and 5+(s) = Caa+(s), (5.2) 

where Ca and are some constants (to be fixed later on) . Notice that from p. 41) it fol- 
lows that the eigenfunctions of f}2(s) are the same functions (|3.2p . but the eigenvalues 
are C^A„ +E, i.e., 

f)2(s)a>„(s) = (c2A„ + £;)$„(s). (5.3) 

A straightforward computation yields 

U^)^b{s)b+is) + ir'-a")Cl + E, (5.4) 
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and 



(5.5) 



+ U^) - Cl{2a{s) + r(.)) + i(| a" t')CI 



or, equivalent ly, 



+ ^cj{s + %){a{s-\) + T{s-\))e^^ 
+ f)i(.)-(2a(.)+r(.))/ + i(§a"-r'). 
The right-hand side of (|5.5p suggests us to use the foUowing new operators 

c{a) = CbKs) e-5^V^(s + 1) = CbCaiais + 1) - e-^^K^)), 
c+(s) = C6v/^(7TT) 65^=6+ (s) = Cfca(a(s + 1) - uis) ), 



(5.6) 



where, as before, i'{s) = \J(j{s + l)[a{s) + t(s)). So, 
[()2(5),c(.)] = -C,2(a"-r')c(s) + CaC, [U^) + {{a"-T')Cl-E) l] a'(s + i), 
[{1,(.),c+(.)] = C,2(a" - r')c+(.) - aaa'(. + i)[f)2(s) + ((^" - r')Q2 - E) /], 
[c(s), c+(.)] = ClCl (a'{s + i)e-^=K«) + ^{s)e^^cj\s + i) - [i.2(.) - _ i)]/) . 
The above expression leads to the following 

Theorem 5.1. If o" — 0, then the operators \)2{s), c(s) and c^(s), defined by 
{5-4^ and h5. respectively, form a closed algebra such that 



[[l2(s), c(s)] = r'C2c(s) + C,Caa'{Q) ([)2(s) - r'C^ - E) , 
[[)2(.),c+(,s)] ^ -r'C2c+(.) - C,Caa'{Q) %{s) - r'C,2 - E) 

[c{s),c+{s)]=Cl[a{s)~a'm^2{s)-E)_. 

Observe also that with this particular choice a'{s + ^) = (t'(0) and 
c(s) + c+{s) = CbCa{^^{s) + 2a'(0) + t{s)), 
v'^is) - v^{s - 1) = (7'(0)(2fT(s) + r(s)) + t'ct(s). 



Furthermore, using the boundary conditions a{s)p{s) 



0, one finds 



h-l 



b-l 



(C$„, $„)<j = CaCb ^ Cr(s + l)$„(s)$„(s) -^V{s~ 1)$„(S - 1)«'™(S) 



s— a 
h-l 



fa-l 



CaCfc ^ a(s + l)$„(s)$™(s) - CaCb ^ J^(s)$„(s)$™(s + 1) 

s— a s — a 

($„,C+$„)d, 
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where 



i.e., the following theorem follows. 

Theorem 5.2. The operators c(.s) and c+(.s) are mutually adjoint. 

Notice also that the operators ()i(s) and i)2{s) arc sclfadjoint operators. 

Remark 5.3. Since A = A„ = -n(T' + (n - l)a"/2), the identity a" = is 
equivalent to the statement that A„ is a linear function of n. In this case A„ = — nr'. 

In the following we will consider only the case when a" = 0, i.e., the case of the 
Meixner, the Kravchuk and the Charlier polynomials. 

If we define the operators 

K_{s) = -t'CI c{s) - CCaa'iQ) {^^{s) - r'Cl - E) , (5.7) 
K+{s) = -r'Cl c+(s) - CtCy{0) ((^^(s) - t'C^ - E) , 

then 

[Ko{s),K±{s)] = ±K^{s) y [K_{s), K+{s)] = AoMs) + A^, 

Ao =-2T'a'{0)C^Ct{-T'Cl){a'{0)+r') and 

A, = -EAoi-r'Cl)-^ + ClCy^a'{0)T{0) - <7(0)t']. 

The case = corresponds to the Charlier case (see the previous section). If 
^0 7^ 0, we have two possibilities: Aq > and Aq < 0. In the following we will choose 
= -l/r', i.e., -t'C'I = 1. 

In the first case Aq > one can choose Cb and E in such a way that Aq = 2 and 
Ai = 0. Thus 

-^/(0)[r'+a'(0)]' 2? " ^^"^^ 

Consequently, the operators K± and Kq are such that 

[Ko{s),K±{s)] = ±K±{s) and [K_{s),K+{s)] = 2Ko{s). (5.9) 

This case corresponds to the Lie algebra Sp(2, DV). 

In the second case one can choose Cf, and E in such a way that Aq = —2 and 
Ai = 0. Thus 

- <7'(0)[r' + a'(0)]' ^ - 2^' " ^^"^"^ 

Consequently, the operators K± and Kq are such that 

[Kois),K±{s)]^±K±{s) and [X+(s), i^_(s)] = 2Xo(s). (5.11) 
This case corresponds to the Lie algebra so(3). 

Notice that since the operator f)2(s) is selfadjoint, the operators K±{s) are mu- 
tually adjoint in both cases, i.e. 
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5.1. Dynamical symmetry algebra Sp(2,JH). Let us consider the first case. 
We start with the operator 

K\s) - Klis) - Ko{s) - K+is)K_is), (5.12) 

where Kq{s), K+{s), and K_{s) are the operators given in (|5.7I) . A straightforward 
calculation gives 

K {s) - E{E - 1)1, E - 2a'(0)(a'(0)+r')' 

where E is given by (|5.8I) . i.e., the K'^{s) is the invariant Casimir operator. 
Furthermore, if we define the normalized functions 



we have 

-l)$„(s), i^o(s)*n(s) = (n + £;)$„(s). (5.13) 

Now using the commutation relation (15. 9p . it is easy to show that 
i^o(s)[i^±(s)$„(s)] = {n + E± l)/^±(s)$„(s). 
Consequently, from (I5.13P and the above equation we deduce that 

K+{s)^n{s) = K„$„+l(s), /C_(s)$„(s) = K„$„_i(s). 

Employing the mutual adjointness of the operators K±, one obtains 

K„ = (ii:+$„(s),$„+i(s))d = ($„(s),is:_$„+i(s))d = K„+i, 

thus 

i4:+(s)$„(s) = K„+i$„+i(s), i^_(s)$„(s) = K„$„_i(s). (5.14) 
In order to compute k„, use (jS.lSp and (|5.14|) : this yields 



E{E -I) = {n + Ef - {n + E)- kI ^ k„ = + 

In this case the functions ($«)" define a basis for the irreducible unitary repre- 
sentation D^[—E) of the Lie group (algebra) Sp(2,9^). 

From the above formula it follows that the functions $„(s) can be obtained re- 
cursively via the application of the operator _R'_|_(s), i.e., 



'^nis) = i i^?(s)<I>o(s), <fo(s) = 

where p(s) is the weight function of the corresponding orthogonal polynomial family 
and do is the norm of the Po(s). 
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5.1.1. Example: The Meixner functions. Let consider the Meixner func- 
tions 



$f (.) = m(^-")/^(1 - /.r/^+^y ^^1^ MZ'^{s), n>0, 

and the hamiltonian [)]^(s) 

f)f (s) = -Vms(s + 7- l)e-^^ - VMs + l)(s + 7)e^= + (s + m (s + 7)) I, 



thus f); (s)$^(s) = (s). In this case we have Cb = E = ^, Ca ^ y 

Therefore 



(g - l + 7)M„-i3, _L + 



s + , 1 



2 ^-^Os 



Consequently, 



Moreover, 



V 1~M V1~A* 



f)2(«) = i)M + 1= b{s)b+is) + 1-1. 



I),(s)$f(s)=(n+J) 
/fo(a) = -V«(s-l + 7)^ e-«» - V(« + !)(« + 7)^ + + 2) /, 



K^is) = - V5iZl±2) - mV(^ + 1)(^ + 7) ^0, + ^ (2, + ^) /, 
1 — /i 1 — M 1 — At 



1 — 1 — M 

and 

K+{s)^^is) = V(n + l)(n + 7)$f+i(s), 

(5.15) 

= ^n(n + 7-l)$^i(s). 
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Using the fact that [)2$o^(s) = |$o^(s), together with the formulas (|^ and (|??T5|) . 
one finds 

= = VK-'^-I + tXCs - 1) - s/.-^<i>o*^(s), 

therefore the normahzed function $o^(s) is 



r(7) Vr(7)r(s + i)' 



and 



n!r(7 + n) 



r(,s + 1) 



A similar result have been obtained before in [8]. 

5.2. Dynamical symmetry algebra so(3). Let us consider the second case 
and define the following operator 



(5.16) 



where Ko{s), K+{s), and K-{s) are the operators given in ()5.7p . Substituting the 
value of E, given by (|5.10p . and doing some straightforward computations yield 



K^{s) ^ E{E -1)1, E^ 



r(0)g'(0) -tV(0) 
2CT'(0)(cr'(0) +r') ' 



i.e., the K'^{s) is the invariant Casimir operator. 

Moreover, if we define the normalized functions as 



dl 



Puis), 



we have 

= E{E - l)$„(s), i^o(s)$«(s) = {n + £;)$„(s). 

Now using the commutation relation (15. lip , we have 

Ko{s)[K±{s)^„{s)] = {n + E± l)A±(s)$„(s). 
Consequently, from (|5.13p and the above equation, we conclude that 

K+{s)<^>n{s) = K„$„+l(s), i^_(s)$„(s) = K„$„_l(s). 

Using the mutual adjointncss of the operators K±, one obtains 
thus 

i4:+(s)$„(s) = K„+i$„+i(s), i^_(s)$„(s) = K„$„_i(s). 



(5.17) 



(5.18) 
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To compute k„, use (I5.13P and (|5.18p : this leads to 



E{E - 1) = {n + Ef + [n + E) + k^^^ k„ = V^(^0(^rr2irnj. 

In this case the functions ($„)„ define a basis for the irreducible unitary repre- 
sentation D'^{—E) of the Lie algebra so(3). 

As in the previous case, from the above formula it follows that the functions <&„(s) 
can be obtained recursively via the application of the operator K^(s), i.e., 



Ki ■ ■ ■ Kn Uq 

where p{s) is the weight function for the associated orthogonal polynomial family and 
do is the norm of the Po{s). 

5.2.1. Example: The Kravchuk functions. Let us consider now the Krav- 
chuk functions 



= p(-")/^(l - pr—^y^^l^^I^K^^is, N), 0<n<N, 



and the corresponding hamiltonian [) ( 



K ^ _ ^/psjN + ^ Np + s-2ps ^ ^p{s + l){N-s) 
"i \^ ~ n ^ ' 1 n ^ 



thus [)f (s)$^'(s) = (s). In this case Ca = ^T"^, ^ Vp^> E - ^ 
therefore 



- -^p{N-s + l)e-2'^^ + v/(l-p)(s + l)e2'^ 



Consequently, 



Moreover, 



Us) - (1 - p) f)i(^) - f - - f - 1- 



ifo(s) = -v/p(l-p)s(iV-s + l)e-^= - ^pil - p){s + 1){N - s)e^ 
+ [7V(p-i)-,s(2p-l)]/, 



ii:+(s) {l-p)y/s{N - s + l)e-^^ +py'is + l){N - 3)6^" - y/p{l-p){2s- N)I, 
K^{s) ^py/s{N - s + l)e-^'' + {l-p)y/{s + l){N - s)e9^ - y/p{l-p){2s - N)I, 
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and 



K^(s)$^(s) = -(iV + 2)<i>«(s), 



N 



i^+(s)<J>f (s) = V(« + l)(iV-n)<I>f+i(s), 

, (5.19) 

i^_(s)$f (s) = ^n(iV-n + l))$^_i(s). 

Using the fact that [)2$^(s) = -f-$^(s), together with the formulas ((5J1) and (|5.19p . 
we find 



^ X_(.)<(.) ^ (^.<(.) - f^\r,^'h ^^{s - l)j , 

therefore the normahzed function $^(s) is equal to 



and 



'^(h) - _ N(jV-«-^)/2 \ n\{N -n)\ 



s / \ 1 — p 



6. The q-case. To conclude this paper we will discuss here briefly what happens 
in the q-case. The preliminary results, related with this case, have been presented 
during the Bexbach Conference 2002 [21- A more detailed exposition of these results 
is under preparation. 

One can first introduce the corresponding normalized functions 



3'n(s) = Pn{s;q), (6.1) 

dn 

where d„ is the norm of the ^-polynomials Pn{s] q), p{s) is the solution of the Pearson- 
type equation 

^ TTWis)p{s)]^T{s)p{s) or a{s + l)p{s + l) ^ cr{-s- fi)p{s), 



and A(s) is an arbitrary continuous function, not vanishing in the interval (a, b) of 
orthogonality of P„. If P„(s; q) possess the discrete orthogonality property (|2.3p . then 
the functions <i>n(s) satisfy 

(<i>„(s),<i>„,(s)) = ^ $„(<,)$„, (s)^^i^ = Sn,m. (6.2) 



Notice that if A{s) = ■\/Vxi(s), then the set ($„)„ is an orthonormal set. Obviously, 
in the case of a continuous orthogonality (as for the Askey- Wilson polynomials) one 
needs to change the sum in (|6.2p by a Riemann integral [SI [55] . 
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Next, we define the g-Hamiltonian ^q{s) of the form 

where 



^,is) ^^A{s)H,is)^, (6.3) 



^cr(-S-/i + l)cr(s) _ ^g(-S-^)o-(s + 1) . 

Vx{s) ^ Ax{s) ^ 

(j{-s-n) cr(s) 



(6.4) 

/\x{s) ' Vx{s) ) ^' 
As in the previous case, one can easily check that 

i3,(s)$„(s) = A„$„(s). (6.5) 
Now we define the a operators: 

Definition 6.1. Let a be a real number and A{s) and B{s) are two arbitrary 
continuous non-vanishing functions. We define a family of a-down and a-up operators 
by 



VVa;i(s) \ 









/ Vx{s) V 


Ax{s) ) 


^A{sy 



(6.6) 



at U) ■= ^ A(s) \ , h^e-^^ - , / ^Li^ ) e"9. 



respectively. The first result in this case is [5]: 

Theorem 6.2. Given a q-Hamiltonian I6'.^|j i^g(s), i/ien i/ie operators Cl^(s) a«d 
aJ^(s), defined in \3.!^) . are such that for all a € C, iOg(s) = ClJ,(s)ci4('S)- 

Our next step is again to find a dynamical symmetry algebra, associated with the 
operator S)q{s), or equivalently, with the corresponding family of g-polynomials. 

Definition 6.3. Let <^ be a complex number, and let a{s) and b(s) be two opera- 
tors. We define the (^-commutator of a and b as 

[a{s),b{s)], = a{s)b{s) - ^6(s)a(s). (6.7) 



We want to know whether the following problem: To find two operators a{s) and 
b{s) and a constant <; such that the Hamiltonian Sjq{s) = b{s)a{s) and [a{s) , b{s)]^ = L, 
has a non-trivial solution. 



Obviously, we already know the answer to the first part: these are the operators 
b{s) = flt (s) and a{s) = Cl'^{s), given in (|6.6I) . The answer to the second part of this 
problem is summarized in the following two theorems (in what follows we assume that 
A{s)^B{s)). 

Theorem 6.4. Let ($„)„ be the eigenf unctions of Sjq{s), corresponding to 
the eigenvalues (A„)„, and suppose that the problem 1 has a solution for A 7^ 0. Then 
the eigenvalues A„ of the difference equation Jg.^p are q-linear or q^^ -linear functions 
of n, i.e., Xn = C'lq" + C3 or A„ = C2(Z^" + C3, respectively. 
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Theorem 6.5. 121 Let S)q{s) he the q-Hamiltonian h6.S\) . The operators h{s) = 
ClJ,(s) and a{s) = Cl;!;,(s), given in i6.6]) with B{s) = A{s), factorize the Hamilto- 
nian S)q{s) i6.S\) and satisfy the commutation relation [a{s) , b{s)]c; = A for a certain 
complex number if and only if the following two conditions hold 



Va;(s) / Va;i(s — l)Vxi(s) a(s — a)a(—s — + a) , 

^ ^, (6. 



Va;i(s — a) y Va::(s — Q;)Aa;(s — a) V (t(s)it(— s — + 1) 
1 / cr(s — a+1) a[~'S — [1, + a)\ 1 / cr(s) (t(— s — /i) 



Aa;(s — a) \^ Va;i(s — a + 1) ' Va:i(s — a) ) ^Vx\{s) \yVa;(s) ' Aa::(s) 

(6.9) 

The proof of the theorem 16.51 is similar to the proof of the theorem 14. 1[ presented 
here for the case of the uniform lattice a;(s) = s. 

Let us point out that the ^(respectively, (7~^)-linearity of the eigenvalues is a 
necessary condition in order to provide that the solution of problem 1 exists. But this 
condition is not sufficient. For example, if we take the discrete g-Laguerre polynomials 
L"(a;; g) (for more details see [2l) with a ^ q^^^^, the problem has not a solution, but 
A„ is a g-linear function of n. 

6.1. Examples. We present here only two examples, others can be found in [5]. 

6.1.1. The Al-Salam & Carlitz I g-polynomials Uri^\x;q). We start with 
the very well known case: the Al-Salam & Carlitz I polynomials j20j. The correspond- 
ing normalized functions p.2[) are 



cfi / ^ , / {qx,a oo -a "g(^) n , / ? 

^n[x) = \ r- — — A^(s)201 

^ q)(q;q)n[q,a,q/a;q)oo \ 




X := q . 



Putting A{s) — •\/Va?i(sy = \/xkq, where kq — q2 — q 2 , we have that the functions 
($„)„ satisfy the orthogonality condition 



f $„(x)$„(a;)|^ 



where the integral f{x)dqX denotes the classical Jackson g-integral. 

For these polynomials a{s) + t(s)Vxi(s) — a, therefore a = and <; = q~^. The 
g-Hamiltonian has the form 



q'^^^a{x — l)(a; — a) _g \J a(l — qx){a — qx) 
[q - 1)'.t2 " ^2 

/ ^{qjx- l)x + a(l + q-qx)) \ _ ^ 

Consequently, i3g(s)$„(s) — gi li-q)i ^n{s) and the operators (a; = q'^] 

1 

aHs) ^ ai(s) = ^ (^^(x-l/q){x-a/q) - l) 
1 

a^s) ^ alis) - f- (V(a;-l)(x-a) e-^= - /) , 



(6.10) 
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are such that 

ai'(s)a-'-(s) and [a^{s),a^{s)]g-i = ^. 

Kg 

A straightforward calculation shows that the operators Cl^(s) and a-''(s) are mutually 
adjoint. A similar factorization was obtained earlier in [5] and more recently in [3] 
with the aid of a different technique. 

The special case of the Al-Salam & Carlitz I polynomials are the discrete q- 
Hermite I hn{x;q), x = q'^ , polynomials, which correspond to the parameter a = — 1 

m- 

6.2. Continuous g-Hermite polynomials. Let us now consider the particular 
case of the Askey- Wilson polynomials when all their parameters are equal to zero, i.e., 
the continuous (/-Hermite polynomials [20]. In this case (t{s) — Caq^^ ■ 

Let us choose A[s) = B{s) = yJVxi{s). In this case a — 1/2 and = 1/q. The 
corresponding Hamiltonian is given by 



1 / Ccq -ds, C,yq 



9),{s) = ...... r^r . e-''+ 



k'^y/sinO \sm{6+^ log g) ^sin(6+i log g) sm{9—^ log g) ^sin(^ — i log q) 



fc| sin 9 \sm{6 + ^ log g) sin(6l — | log g) 
and the a-operators 



r< „2s 1 / r' n-^s 



= e2"% . „ . „ --e 2 



fc^ sin 9 sm{6 + | log q) V ^'^^ ^ sm{6 — ^ log q) ' 



"■i/^^"^' y -fc2 sine* sin(0 + flog g)"" Y -/c2sin6lsin(6l- f logg)^ ' 

are such that 

4C 

a^{s)a^{s) = Sjq{s) and [a^{s),a'^{s)]i/q = 

Kq 

Another possible choice is A{s) = B{s) = 1 [TT|, hence a straightforward calcula- 
tion shows that the two conditions in Theorem 16.51 hold if <^ = g^"'^, thus A — ACcrk^^. 

With this choice the orthogonality of the functions is /^^ $„(s)3>m(s)da; = Sn,m- 
In this case, the Hamiltonian is 



fc2 I sinesin(6'+f Ing) sin(6l- 1 In g) sin 6* g + g"^ - 2 cos 261 

and 



(6.11) 
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With these operators 

4C 

S)q{s) = a^{s)a^{s) and [a^{s),a'^{s)]q-i = — 

ftq 

This case was first considered in [TT]. 

Acknowledgments. We are grateful to A. Ruffing for inviting us to participate 
at the Bexbach Conference 2002 and present there our results The research of 
RAN has been partially supported by the Ministerio de Ciencias y Tecnologfa of 
Spain under the grant BFM-2003-06335-C01, the Junta de Andaluci'a under grant 
FQM-262. The participation of NMA in this work has been supported in part by the 
UNAM-DGAPA project INI 12300. 

REFERENCES 

[1] R. Alvarcz-Nodarsc, Polinomios hipergemetricos y q-polinomios. Monograffas del Seminario 
Matematico "Garcia de Galdeano" Vol. 26. Prensas Universitarias de Zaragoza, Zaragoza, 
Spain, 2003. 

[2] R. Alvarez-Nodarse, N. M. Atakisliiyev, and R. S. Costas-Santos, Factorization of 
hypergeometic-type difference equations on nonuniform lattices: algebraic properties of 
the raising and lowering operators (in preparation). 

[3] R. Alvarez-Nodarse and R,.S. Costas-Santos, Factorization method for difference equations of 
hypergeometric-type on nonuniform lattices. J. Phys. A: Math, Gen. 34 (2001), 5551-5569. 

[4] R. Askey and S. K. Suslov, The g-harmonic oscillator and an analogue of the Charlier polyno- 
mials. J. Phys. A. 26 (1993), L693-L698. 

[5] R. Askey and S. K. Suslov, The g-harmonic oscillator and the Al-Salam and Carlitz polynomials. 
Lett. Math. Phys. 29 (1993), 123-132. 

[6] N. M. Atakishiyev, Construction of the dynamical symmetry group of the relativistic har- 
monic oscillator by the Infeld-HuU factorization method. In Group Theoretical Methods in 
Physics, M. Serdaroglu and E. Inonu. (Eds.). Lecture Notes in Physics, Springer- Verlag, 
180, 1983, 393-396; Ibid. Theor. Math. Phys. 56 (1984), 563-572. 

[7] N. M. Atakishiyev, A. Frank, and K. B. Wolf, A simple difference realization of the Heisenberg 
g-algebra. J. Math. Phys. 35 (1994), 3253-3260. 

[8] N. M. Atakishiyev, E.I. Jafarov, S.M. Nagiyev, and K.B. Wolf, Meixncr oscillators. Revista 
Mexicana de Fisica 44(3) (1998), 235-244. 

[9] N. M. Atakishiyev, M. Rahman, and S. K. Suslov, On Classical Orthogonal Polynomials. Gon- 
structive Approximation. 11 (1995), 181-226. 
[10] N. M. Atakishiyev and S. K. Suslov, Difference analogs of the harmonic oscillator. Theor. Math. 

Phys. 85 (1991), 442-444. 
[11] N. M. Atakishiyev and S. K. Suslov, A realization of the (/-harmonic oscillator. Theor. Math. 

Phys. 87 (1991), 1055-1062. 
[12] N. M. Atakishiyev and B. Wolf, Approximation on a finite set of points through Kravchuk 

functions. Revista Mexicana de Fisica 40 (1994), 1055-1062. 
[13] G. Bangerezako, Discrete Darboux transformation for discrete polynomials of hypergeometric 

type. J. Phys A: Math. Gen. 31 (1998), 2191-2196. 
[14] G. Bangerezako, The factorization method for the Askey- Wilson polynomials. J. Gomput. Appl. 

Math. 107 (1999), 219-232. 
[15] C. Berg and A. Ruffing, Generalized g-Hermite polynomials. Gomm. Math. Phys. 223 
(2001), 29-46. 

[16] L. C. Biedenharn, The quantum group SUq{2) and a q-analogue of the boson operators. J. 

Phys. A. 22 (1989), L873-L878. 
[17] R. Floreanini, J. LeTourneux, and L. Vinet, More on the g-oscillator algebra and g-orthogonal 

polynomials. J. Phys. A: Math. Gen. 28 (1995), L287-L293. 
[18] R. Floreanini and L. Vinet, g-orthogonal polynomials and the oscillator quantum group. Lett. 

Math. Phys. 22 (1991), 45-54. 
[19] L. Infeld and T. E. Hull, The factorization method. Review of Modern Physics 23 (1951), 21-68. 
[20] R. Koekoek and R. F. Swarttouw, The Askey-scheme of hypergeometric orthogonal polynomials 

and its q-analoguo. Reports of the Faculty of Technical Mathematics and Informatics No. 

98-17. Delft University of Technology, Delft, 1998. 



18 FACTORIZATION OF DIFFERENCE EQUATIONS 



[21] A. Lorek and A. Ruffing, and J. A. Wess, A qr-deformation of the harmonic oscillator. Z. Phys. 
C 74 (1997), 369-377. 

[22] M. Lorcntc, Raising and lowering operators, factorization method and differential/difference 

operators of hypcrgcomctric type. J. Phys. A: Math. Gen. 34 (2001) 569-588. 
[23] W. Miller (Jr.), Lie theory and difference equations. 1. J. Math. Anal. Appl. 28 (1969), 383-399. 
[24] W. Miller (Jr.), Lie theory and g-differcnce equations. SIAM .J. Math. Anal. 1 (1970), 171-188. 
[25] Sh. M. Nagiyev, Difference Sehrodinger equation and g-oseillator model, Theor. Math. Phys. 

102 (1995), 180-187. 

[26] A. F. Nikiforov, S. K. Suslov, and V. B. Uvarov, Classical Orthogonal Polynomials of a 
Discrete Variable. Springer Series in Computational Physics. Springer- Verlag, Berlin, 1991. 

[27] A. F. Nikiforov and V.B. Uvarov, The Special Functions of Mathematical Physics, Birkhauser, 
Basel, 1988. 

[28] Yu. F. Smirnov, Factorization method: New aspects. Revista Mexicana de Fisica 45 (1999), 1-6. 



